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ABSTRACT 

The collapse of marginally Jeans unstable primordial gas clouds in the presence of 
UV radiation field is discussed. Assuming that the dynamical collapse proceeds ap- 
proximately in an isothermal self-similar fashion, we investigate the thermal evolution 
of collapsing central core until H2 cooling dominates photoheating and the tempera- 
ture drops to below lO^K. Consequently, the mass of the cooled core is evaluated as 
Mcoo\ = 3.6 X 10^ Mq (/2i/1)~°''^^- This scale depends only on the incident UV inten- 
sity, and provides a lower limit to the mass of collapsed objects in the UV radiation 
field. 
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1 INTRODUCTION 

Formation of primordial objects, such as young galaxies and 
globular clusters, in the early universe is a fundamental 
problem in modern cosmology. Rapid progress of observa- 
tions in the last decade has pressed theoreticians to con- 
struct the physically correct and proper theory to solve this 
problem. Recently, formation of primordial objects has been 
investigated mainly in the following two contexts; one is the 
formation of "first stars" or "first luminous objects" (e.g. 



( Pegmark et al. 1997 ; i 


Vbel et al. 1995 ; 


Omukai&Nishi 1998 


pSfishi & Susa 19991) ), and the other is that of "second gen- 


eration objects" (e.g. 


Haiman, Rees & Loeb 1997 


Kepner. 



Babul fc Spergel 1997] pmukai fc Nishi 19991 ) ). These two 
populations are expected to arise in quite different physical 
environments. For instance, the former is likely to be born 



from the 



pimiuidial gas uuUei ULLle mimeuces uf Llie ex Lei - 
nal rad iation field, except for that of the ooamic microwave 



backgrtiuud ladiaLiuu (CMD). The laLLe 



jrcT 
get 



Llie uLlitii liauu 



is largel y nffGctcd by the external UV radiation produced by 
the former. The UV photons not only ionize and heat up the 
pregalactic gas clouds, but also dissociate H2 in the clouds, 
which is an important coolant in metal-free gas clouds. In 
this paper, we pay particular attention to the formation of 
the latter population. 

At high redshifts, say z S> 5, ionized bubbles around 
photon sources, such as young galaxies or AGN's, are still 
so small that most of the gas in the universe is not ion- 
ized and photoheated. In this case, the external radiation 
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affects the form ation of primordial objects q nlv via H2 pho- 
todissociation (Haiman, Rees & Loeb 1997). On the other 
hand, gas clouds near the ionizing sources, or those at lower 
redshifts {z ^ 5) are ionized and photoheated, before they 
collapse and cool. Therefore, the formation of population 
III objects close to the ionizing sources and the formation 
of galaxies at low redshifts are likely to start from hot ion- 
ized media. Kepner, Babul & Spergel (1997) investigated 
this problem in the context of dwarf galaxy formation un- 
der the assumption that pregalactic clouds are in hydro- 
static equilibrium. They found that the gas transit from H"^ 
phase to H and H2 ones, as radiative cooling proceeds. Cor- 
belli, Galli & Palla (1997) also found similar phenomena in 
the rotationally supported hydrostatic gaseous disc. How- 
ever, gravitat ional collapse of ionized gas proceeds almost 
isoth e rmally (Umemura fc Ikeuchi 1984 iThoul fc WeinbergI 



199^ ; [Kitayama fc Ikeuchi 2000| |Susa fc Umemura 2000|) 



and such a collapse becomes inevitably dynamical (Larson 



196S| ). 

In this paper, we investigate H2 cooling in the dynam- 
ically collapsing core in marginally Jeans unstable clouds. 
If these clouds are self-gravitating, they are likely to col- 
lapse almost spherically with the temperature kept nearly at 
~ lO^K, until the external radiation field is shielded by the 
clouds themselves. We thus emplov the i sothermal Larson- 
Penston similarity solution (Larson 1969) and solve explic- 
itly non-equilibrium chemical reactions and energy equation 
to trace the thermal history of a collapsing central core. 
Consequently, the mass and the size of the cooled core are 
assessed, and their astrophysical implications are discussed. 
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2 H2 COOLING IN THE COLLAPSING CORE 



2.1 Assumptions and Equations 



present case and the spectral change in the Lyman- Werner 
bands will not alter our results |^ 

The photoionization rate is computed taking account 



• 1 . 1 11 r -ir i- • 1-1 of the absorption of ionizing photons as ( Pajiri & Umemura 
We con nidor tho oollapao ot a aolt gravitatmg primordial gao 1 1 — ^1 — i -J 



cloud e xposed to the isotropic UV radiation field. If the UV 



1998 



intensity is strong enough to penetrate through the cloud, 
the temperature of the cloud is first kept at ~ 10*K, owing 
to the balance between photoheating and radiative cooling. 



in a self-similar fashion ((Jmemura & Ikeuchi 1984; Thoul & 


Weinberg 1996 




Kitayama & Ikeuchi 200(\, Susa & Umemura 



200c). The isothermal collapse of a self-gravitating system 



always te nds to conve rge to the Larson-Penston similarity 



solution (Larson 1969) because it is a relatively stable so- 



lution (e.g. dHanawa fc Nakayama 1997 )). We thus assume 



that the cloud obeys the Larson-Penston solution, which is 
characterized by the central density 



Pc{t) = 



Pc(0) 



(t Alp (0) 



(1) 



and the core radius 



rdt) = (tLP (0) -t)c. 



(2) 



Here t denotes the elapsed time from an initial instant char- 
acterized by pc(0). Ca is the sound speed in the core, and t^p 
is defined as tLP (0) = -\/l.667/47rGpc (0) In order to elu- 
cidate the thermal history of nearly isothermally collapsing 
core, we employ an one-zone approximation and integrate 
numerically the following energy equation: 



dec 
'dt 



Pcdpc 
p§ dt 



Pc 



(3) 



where Pc is the central pressure, ec is the specific energy per 
unit mass, and A and F denote the cooling and heating rates, 
per unit volume, respectively. The cooling rate includes H2 
cooling by rovibrational transition lines, as well as H atomic 
cooling. We also solve time-dependent non-equilibrium re- 
actions for e, H, H+, Hj, H~, and H2, to evaluate A. Unless 
stated explicitly, reaction rates are taken from the recent 
compilation by Galli & Palla (1998). H2 photodissociation 
rate is evaluated using the self-shielding function in Draine 
& Bertordi (1996). photodissociation is also taken into 
account, based upon the cross section in Standi (1994). H~ 
radiative detachment is assessed using the fitting formula 
for the cross section in Tegmark et al.(1997). 

The incident UV intensity is set to be PJ^ = 
/2i(i'/^'i)~^10~^^erg s~^cm~'^str~^Hz~^, where vl is the 
frequency at the Lyman limit. We take various values for 72i 
in the range 10~^ hi ^ 10^. The power-law approxima- 
tion of the UV spectrum adopted here is admittedly rather 
crude. In reality, the UV intensity in the Lyman- Werner 
bands could be modulated by so-called the "sawto othing ef 



feet" and reduced by a factor of 10 — 10 at 2: = 15 ( Haiman 



Abel & Rees 2000). As will be discussed later, however, the 
two-step photodissociation process is not important in the 



Susa fc Umemura 2000) 



-r—(Jv dQ.dv 
hu 



h 



3r.,(n)4/3 



7(4/3, r,^ {n))dn 



47r/i" JiKKT^i 7(4/3,T^^) 



4/3 



(4) 



where (t^ is the photoionization cross section which is pro- 
portional to ly'"^, O, is the solid angle, and 'y{a, b) represents 
the incomplete gamma function. The optical depth at the 
Lyman limit is defined as r^^^ = nmrcO'i.j^, where rc is as- 
sessed by equation (^) and nm denotes the central HI num- 
ber density. The accreting envelope outside Vc only alters 
the optical depth t^^ less than a factor of 2. Similarly, the 
photoheating rate is expressed as 



r = nm 



hu 



twL)dQ,du 



47r in 
3 



7(l,r.J 7(4/3,r.J 



erg s cm 



(5) 



The initial central density is taken as twice the self- 
shielding critical density ricr derived in Tajiri & Umemura 
(1997), where ricr represents the density above which the 
HI fraction exceeds 0.1. Thus, the numerical integrations 
of the energy equation and chemical reactions start when 
the cloud center begins to be almost neutral. The initial 
temperature and chemical abundances are fixed assuming 
thermal and chemical equilibria. Even if they are perturbed, 
they immediately converge to the equilibrium values. We 
stop the numerical integrations when the temperature of 
the core drops to below 5000K and the assumption of an 
isothermal cloud breaks down. 



2.2 Numerical Results 

The numerical results are presented in Figs|l| and^. In Fig.^, 
time evolution of the core temperature (upper panel) and 
the chemical compositions (lower panel) are plotted for four 
different values of the incident intensity. In every case, the 
temperature evolves nearly isothermally after the collapse 
calculation starts. This ensures that the density evolution 
approximately follows the isothermal Larson-Penston simi- 
larity solution. After the isothermal collapse, core temper- 
ature drops dramatically, with the time-scale much faster 
than the collapse time, and the loci on the pc — Tc plane be- 
come nearly vertical. This corresponds to the transition from 

■t- The sawtoothing effect could also suppress destruction of by 
the UV photons, which could in turn enhance the H2 abundance. 
The assumption of a power-law spectrum should thus provide the 
limiting case in which H2 formation is maximally suppressed. In 
practice, however, this has only minor impacts on our present 
results (see also Section 3). 
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the H phase to the H2 one mentioned in Kepner, Babul and 
Spergel (1997). The increased cooling rate by H2 overwhelms 
the photoheating rate which is reduced due to strong self- 
shielding. In the lower panel of Fig.|l| fractions of electron 
and H2 are plotted. At T ^ 10* K, the electron fraction is 
kept at ~ 10"^ by ionizing UV photons. At T ^ lO"* K, the 
photons are self-shielded and recombination proceeds. How- 
ever, the radiative cooling time becomes shorter than the 
recombination time, and the electrons become out of equilib- 
rium. Significant amount of free electrons thus "freeze out" 
as the system cools via Ly-a and H2 ro-vibrational tran- 
sitions. The ionization degree still remains at the level of 
~ lO"'^, even at T ~ 5000K. Feeded these relic electrons, 
H2 molecules are formed to the level of ~ 10~^. Remark 
that the similar phenomena h ave been found in the post- 
shock r egion of p rimordial gas ( Shapiro fc Kang 1987 ; Kang 



1998 



Nishi et al. 1998) 



& Shap iro 1992 ; Susa et al. 

In Fig.^ various time-scales are plotted for I21 ~ 0.1 
and 10^ The upper panel shows time-scales related to the 
energy equation, i.e., the collapse time (tLp), the cooling 
time (tcooi), the II2 cooling time (tcooi,H2)' ^^'^ pho- 
toheating time (tuv). It is clear that temperature is de- 
termined by thermal equilibrium, because the condition 
tuv ~ fcooi is well satisfied. The transition from the H phase 
to the H2 one takes place as H2 cooling dominates UV heat- 
ing. We remark that the collapse time-scale ti,p is much 
longer than the other thermal time-scales. Therefore, hy- 
drostatic equilibrium assumption is a poor one owing to the 
cooling instability. 

In the lower panel of Fig.^ time-scales related to H2 
formation/destruction are plotted. Chemical equilibrium is 
achieved throughout the evolution because the H2 dissocia- 
tion time is almost equal to the H2 formation time {tdis — 
tfor). It is clear that two-step photodissociation (Solomon 
process; isoi) is not important during this calculation. Dis- 
sociation of II2 is dominated by a coUisional process ( II2 + 
at such high temper ature (|Corbelli, GaUi 



H+ ^ H+ + H 



Falls 1997; gusa fc Umemura 200C ) 



3 PROPERTIES OF COOLED CORE 

As shown in the previous section, the runaway collapsing 
core cools rapidly as soon as H2 cooling becomes effective. 
The core will probably cool down to ~ 100 K, which is the 
lowest temperature achievable by H2 cooling. The cooled 
core will then be a free-falling sphere, since the cooling time- 



scale is much shorter than the collapse time (Fig.2). Con- 
sequently, violent collapse ot tlie cooled core will lead to 
subsequent star formation in the cloud center. 

In Figj^, the mass and the radius of the cooled core 
are plotted against the incident UV intensity. The upper 
panel shows the cooled mass, which is approximately fit- 
ted by Mcooi = 3.6 X 10*^ Mq (/2i/1)"'''^^ Remark that this 
scale depends only on the incident intensity /21. This de- 
pendence can be understood essentially as follows. Mcooi is 
roughly estimated from three conditions at the onset of H2 
cooling: 1) Mcooi ^ M](T = 10*K), 2) F ~ Au^, and 3) 
yu2 — 10"* — 10"'^. Condition 1) simply dictates that the 
core mass is comparable to the Jeans mass. Condition 2) 



the II2 fraction (denoted as yH2) at this moment roughly 
converges to a fixed level (Fig.nl) H. These three conditions 

— 1/3 

yield Mcooi ex. I21 , which is slightly steeper than our nu- 
merical results. This difference is caused by two reasons. The 
first one is the destruction of and H" by the UV flux 
which are neglected in the above simplified estimation. For 
large J21, the fractions of and H" are reduced, and the 
amount of II2 produced and II2 cooling rate both become 
smaller. In order to achieve the balance between cooling rate 
and photoheating rate, the system has to collapse even fur- 
ther. As a result, the cooled mass becomes slightly smaller 
than the analytic estimation and the Mcooi — /21 correlation 
becomes steeper for large 72i. The other reason is the rel- 
atively short free-fall time compared to the cooling time at 
lower density. For low 72i, H2 cooling starts to operate at 
low density (Fig.|l]) when the difference between the free-fall 
time and the cooling time scale is small. This makes gentler 
the drop of the temperature as a function of density at the 
onset of H2 cooling (Fig.^, upper panel) . Consequently, the 
cooled mass becomes larger and the A/cooi — ^21 correlation 
becomes shallower for small /21. 

In the limit of very strong UV field, II2 formation is 
completely suppressed as H" and are both destroyed by 
the UV photons. We find that this happens at the threshold 
intensity of /21 — 3.6 x 10* in the case of a = 1. 

The radius of the cooled core is also plotted for various 
721 in the lower panel. It also has the same dependence on 
721 as Mcooi, '"cool tx 7^^*^ "^^. Note that the absolute values 
of cooling radius shown in Fig.^ are larger than the Jeans 
length of the cooled objects, since they are assessed just at 
the onset of II2 cooling. In reality, the temperature drops 
rapidly from ~ 10*K to ~ lOOK due to H2 cooling. There- 
fore, the corresponding Jeans length should be reduced by 
a factor ~ 10^, provided the cooled core mass is constant 
during the free-fall collapse. This Jeans length will provide 
the actual size of the cooled object. 

For 1 ^ 721 ^ 10^, the cooled mass is ~ lO'^ - IO^'Mq. 
Such high intensity is realized near the luminous Population 
III objects, such as young galaxies or quasars. For instance, 
721 = lO'^ roughly corresponds to the intensity at lOkpc 
away from the center of active galactic nuclei whose lumi- 
nosity is 10** erg s~^. In this case, the cooled mass in Fig.^ 
should provide a lower limit to the mass of the cooled ob- 
jects around or inside the young galaxies/quasars. The mass 
(~ 10''A7q), compactness (~ a few pc, 100 times smaller 
than rcooi in Fig.^) and location (~ 10 kpc) of the cooled 
core roughly acc ount for those of globular clusters ( Binney & 



Merrifield 1998). This result implies that the globular clus- 



ters might be formed in the runaway collapsing core around 
luminous host galaxies. 

For 10"^ ^ 721 ^ 1, the cooled mass Mcooi is ^ 



10' 



10 M0, and cooling radius rcooi ^ Ikpc. These in- 



means that the thermal equilibrium is achieved until the 
onset of II2 cooling. Condition 3) is based on the fact that 



3 J/H2 ^t the onset of H2 cooling is basically determined by chem- 
ical equilibrium with given temperature and electron fraction. 
As the temperature drops to T ~ 8000K, the dominant cooling 
mechanism switches from H Ly-Q cooling to H2 cooling. At this 
temperature, the electron fra ction converges to Wrl^^MZl 

re- 

gardless of initial conditio ns ( Susa ct al. 199^ ; Kitayama, Susa, 



Umemura fc Ikcuchi 200o| ) and the H2 fraction freezes out nearly 



at its equilibrium value. 
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tensities correspond to the UV background radiation field 
i nferred from so-called the proxim i ty effect of Lya forest s 
( iBajtlik, Duncan fc Ostriker 1988| ; piallongo et al. 1996| ). 



The actual mass of the cooled gas is determined by some 
mechanism (SN, radiation feedback) which halts the mass 
accretion from the envelope onto the cooled core. In this 
sense, these mass scales provide a lower limit to the mass 
of coo led objects under the UV background radiation field . 
HST (|Pascarelle et al. 1996|) and SUBARU ([Yamada 199^) 



in fact found the "building blocks" which are very compact 
Ikpc), although the mass is still unsettled. The creation 
mechanism of low mass compact objects discussed in this 
paper might be able to explain the formation of such ob- 
jects. Finally, let us remark on influences of dissipationless 
dark matter on the present results especially in case we ap- 
ply them to the formation of dwarf galaxies. If the dark 
matter dominates the gravity of the collapsing core, the col- 
lapse will not follow the Larson-Penston similarity solution. 
The dark matter distribution at the center of a pregalactic 
cloud is still highly uncertain, but the cooled mass in Figj^ 
is likely to be reduced by the dark matter gravity. This point 
will be discussed in more detail in our futu re publications 
(Kitayama, Susa, Umemura & Ikeuchi 200C). 
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Figure 1. Time evolution of tlie physical quantities in the collapsing core for four different values of the incident UV intensity 721- The 
upper panel shows the evolution of temperature, while the lower panel that of electron fraction (solid lines) and H2 fra<;tion (dotted 
lines). 
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Figure 2. Various time-scales during the evolution for /21 = 0.1 and 10^. The upper panel shows the collapse time of the Larson-Penston 
similarity solution (solid line; tLp), the total cooling time (short dased; tcooOi the H2 radiative cooling time (long dashed; tcool,H2)' 
and the photoheating time (dotted; tuv)- the lower panel, time-scales related to H2 formation and dissociation are plotted; the H2 
photodissociation time (dotted; tsoOi the H2 coUisional dissociation time (solid; tdis)i and the H2 formation time (short dashed; tfor)- 
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Figure 3. Mass (upper panel) and radius (lower panel) of the core at Tcore = 5000 K, i.e. at the onset of H2 cooling, as a function of 
the incident UV intensity l2i- 
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